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ABSTRACT 


A code to simulate material transport through porous media was developed at 
Oak Ridge National Laboratory. This code has been modified and adapted for 
use at Lawrence Livermore Laboratory. This manual, in conjunction with report 
ORNL-4928, explains the input, output, and execution of the code on the 
Octopus Computer Network. 


INTRODUCTION 


A transport code developed by Duguid and Reeves! has heen adapted to the 
Octopus Computer Network. The cade calculates two-dimensional (vertical 
section) transport of a dissolved constituent hy ground water. The code has 
been used to calculate transport in regional ground-water systems,” These 
studies were carried out (under contract to the Nuclear Regulatory Commission) 
to investigate geohydrological criteria important in the geological disposal 
of nuclear waste. However, the code can be applied to other ground-water 
transport problems. This manual is designed to quide the user in the input 
and output operations of the code and in the execution procedures required for 
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use at a remote terminal. Duguid and Reeves’ is wel? documented and should 


he used in conjunction with this manual. 


This code is designed to calculate the concentrations of a material within 
porous media. Among its capabilities, it considers transient or steady-state 
flow, saturated or unsaturated conditions, advective transport, hydrodynamic 
dispersion, and radioactive decay. The theoretical development of the 
governing equations and numerical implementation are found in Duguid and 


Reeves, 


en A 


DATA FILE PREPARATION 


Construction of the finite element mesh and generation of a ground-water flow 
field 3,4 precedes preparation of the data file, but the flexibility in this 
procedure puts it beyond the scope of this manual. What follows is a list of 
data requirements. This list should be used in conjunction with ORNL-4928. | 


(7) PROBLEM NUMBER AND IDENTIFICATION LINE 


VARIABLE NAME IDENTIFICATION FORMAT (15,9A8) 
NPROB - Number of the problem 
TITLE - Title of the problem 


(2) CONTROL LINE 


VARIABLE NAME IDENTIFICATION FORMAT( 1416) 
NNP - Number of nodal points 
NEL - Number of eiements 
NMAT - Number of materials 
CM - Number of elements with corrected material properties 
NTI - Number of time increments (iterations) 
NBC - Number of constant Dirichlet boundary conditions 
NST - Number of element sides on which flux-type Neumann 
boundary conditions are applied. 
NTST - Number of elements having seepage surfaces 
KVI - Darcy velocity input control 
KSTR - Control parameter for storage of output 
KSS - Steady state control 


ISSENT - If ISSENT=0, TSFXXXXX is calculated. If ISSENT=1, 
TSFXXXXX exists from a preceding calculation and js 
automatically read. 


(3) TIME CONTROL LINE 


VARIABLE NAME IDENTIFICATION FORMAT(15,5F 10.0) 

IPLOT - Control parameter for saving concentration for 
postprocessing 

DELT - Time increment 

CHNG - Multiplier for increasing the time increment 

DELMAX - Maximum value of DELT 

TMAX - Maximum simulation time 

W - Time-integration parameter 


If [PLOT<0, the code will save concentration values at 20 equally spaced 
values of TMAX. If IPLOT=0, the user must provide the times where the code is 
to save concentrations, The maximum number of times entered cannot exceed 

20. These numbers are read with 5F10.0 format. If IPLOT>0, concentration 
values will be saved at TIME=IPLOT and at TIME=TIME+JPLOT with a maximum of 20 
sets of time saved. 


(4) OUTPUT CONTROL PARAMETER LINE 


VARIABLE NAME IDENTIFICATION FORMAT({I1) 
KPRO - Printer control for steady state and initial conditions 


(5) SOURCE NODE LINE 


VARIABLE NAME IDENTIFICATION FORMAT( 1) 
II - Number of source nodes that can be turned off during 
simulation, 


If the value of II is zero, no source nodes will be turned off and the source 
node contro] line (next input line) should not be entered. Also, the code 
will adjust all source nodes or none at all. 


Aah 


(6) SOURCE NODE CONTROL LINE 


VARIABLE NAME IDENTIFICATION FORMAT(2F 10.0) 
TIMESOURCE - This variable is the simulation time when source 
nodes can be turned off or diminished. 
SOURCE - This variable is the value (strength) of the source 


node after it is adjusted to the new source value. 


In the most general case, the number of cards needed here is equal to the 
number of source nodes. No cards are needed if none of the source nodes are 


adjusted. 
{7) MATERIAL PROPERTIES LINES 


In the most general case, NMAT sets of parameters are required, 


VARIABLE NAME IDENTIFICATION FORMAT(8F 10.0) - first 

line 
FORMAT(F10.0) - second 
line 

KD - Distribution coefficient of an element 

RHOB - Bulk density 

AL - Longitudinal dispersivity of material 

AT - Transverse dispersivity of material 

LAMBDA - Radioactive decay constant 

POR - Porosity of an element 

ALP - Modified coefficient of compressibility 

AM - Molecular diffusivity 

TAU - Tortuosity 


The above material properties values are read into the array PROP, then stored 
as the variable names listed above. Because of this they must be entered in 
floating point form. 


(8) LINES DEFINING NODAL POSITION AND DIRICHLET CONCENTRATION-TYPE BOUNDARY 


CONDITIONS 


In the most general case, NNP lines are required, 


VARIABLE NAME 


IDENTIFICATION FORMAT(15,12,13,3F 10.0) 


M(temporary) - Number of the nodal points 


IM - 


NPN - 


Ys 


Z - 
BB - 


Variable for automatic generation of Dirichlet 
concentration boundary 

Contains node number where constant Oirichlet boundary 
conditions occur 

X coordinate of node 

Z coordinate of node 

Constant material concentration for Dirichlet boundary 
condition 


(9) ELEMENT DEFINITION LINES 


VARTABLE NAME 


IDENTIFICATION ye FORMAT(6[5) 


M(temporary) - Element number 


TE(1,M) - 
JE(2,M) - 
1E(3,M) 
IE(4,M) 
1E(5,M) 


Node I 
Node J 
Node K 
Node L 
Material identification 


NEL lines are usually required. Element cards must be in element sequence. 
If element cards are omitted, the program automatically generates the omitted 


information by increasing the node values of the preceding element by one. 
The material identification is the same as the preceding elemenc. The last 
element card must always be supplied. Nodal sequencing I,J,K,L is 
counterclockwise around the element, starting at Tower left corner of element. 


(10) MATERIAL CORRECTION LINES 


VARIABLE NAME IDENTIFICATION FORMAT ( 1416) 
MI - Element number 
MTYe - Material type 
MK - Element number 
MINC - Increment in element number 


(11) INITIAL CONDITIONS LINES 


In the most general case there is one card per node, i.e., a total of NNP 


cards. 
VARIABLE NAME IDENTIFICATION FO?4AT(15,5X,F 10.0) 
NJ - nodal? point number 
P>(NJ) - Pressure head at current time step 


If groups of adjacent nodal points NJ have identical values RP(NJ), a data 
omission can be used for automatic generatiun of these values. The initial 
concentration of each succeedi1g node will have the same concentration value 
until a new value is encountered. 


(12) NEUMANN FLUX-TYPE BOUNDARY CONDITIONS LINES 


VARIABLE NAME IDENTIFICATION FORMAT (315,5X,2F10.0} 
NI - Nodal point number 
NJ - Nodal point number 
KINC - Increment contro! used for automatic generation of 
boundary conditions 
EJ - Flux at NI normal to the surface (NI,Nd) 
Ed - Flux at NJ normal to the surface (NI,NJ) 


If KINC>0, then the nodal-point increment is ford r7om NI and NJ of the 
immediately preceding card: 


NPINC = Nd - NI 

Two sequences are formed: 
NI + NPINC, NI + 2*NPINC... 
NJ + NPINC, NJ + 2*NPINC... 


Both are terminated when the largest integer is reached that is less than both 
current values of NI and NJ. Corresponding nodal points for these two 
sequences define a surface. Quantity EI is the dot product of the flux at NI 
with an outwardly directed unit vector normal to the element side (NI,NJ). A 
similar definition holds for Ev.) 


(13) SEEPAGE SURFACE ELEMENT LINES 


VARTABLE NAME IDENTIFICATION FORMAT (1416) 
NTSE({MP) - Element having seepage surfaces 
IS(MP,1) - Surface identification array for seepage elements 
IS(MP,2) - Surface ider:tification array for seepage elements 
KINC - Increment control used for aitomatic generation of 
boundary conditions. If KINC is greater than 0, 
automatic generation is employed. 


(14) INPUT LINES OF DARCY VELOCITY AT TIME = Q, z 


The Darcy velocities may be entered in the user's data deck or by reading the 
binary format file VELOUT generated by the TAYLOR ground-water code (GWFIOW). 


The method used to enter the velocities is determined by the vatue of variable 
KVI entered on line (2) of data deck. 


If KVI<O, NEL velocity cards are required 

If KVI=0, no Darcy velocities are read 

If KVI>0, Darcy velocities are read from file VELOUT generated by 
TAYLOR fiow code (GWFLOW). 


The Darcy velocity data deck is read in the following form: JQ is an index 
from 1 to 4 and M is the element number: (VX(JQ,M),JQ=1,4, VZ(JQ,M),JQ=?,4). 


WX(1,M) Velocity in X-direction of node I, element M 
VX(2,M) Volocity in X-direction of node J, element M 
VX(3,M) Velocity in X-direction of node K, element M 
VX(4,M) Velocity in X-direction of node L, element M 
"Z(1,M) Velocity in Z-direction of nede I, element M 
vz(2,M) Velocity in Z-direction of node J, element M 
VZ(3,M) Velocity in Z-direction of node K, element M 
¥Z(4,M) Velocity in Z-direction of node L, element M 


Wodes I,J,K,L run counterclockwise starting with node I at the lower left 
corner of the element. If velocity values are read from cards, format(8E12.4) 


is used. 


PRESSURE HEAD {H{I)) and TOTAL HEAD (HT{1}) at TIME=0 are set to zero 
internally by the code. If the type of problem desired to be run dictates 
that the TOTAL HEAD and PRESSURE HEAD have some value other than zero, the 
code will have to be modified accordingly. The two dimensional array TH(1,1) 
to TH(4,NEL) contains the water content. At present, the water content is set 
(by the code) equal to the porosity of the materiat. 
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RUNNING THE TRANSPORT CODE FROM TTY 


The transport code is a public file (MATTRAN) on the CDC 7600 computers of the 
LTSS computer network. The following describes the teletype interactions. 


USER: MATTRAN / t v 

TTY: BOX NUMBER AND INPUT FILE NAME 

USER: AXX FILENAME Format(A3, 1X,A10) 
TTY: PROBLEM INFORMATION INPUT 

USER: XXXXK J Format(AS,1X,A1) 


The code now echos the user's input. 


ENPUT FLLE se: sac “a iad Sac Ya ea eat FILENAME 
BOX NUMBER... 2... 2.20. AXX 
RESTART DUMP FILE ....-. XXXXX 
HSP SENTINEL... 2 ee eee J 


If all is o.k., type yes; if not, hit return and retype. 
USER : YES (or carriage-return) 


TTY : ALi DONE 


In the above example, the terms used are defined as follows: 
(1) MATTRAN - The transport controllee; a public file 


(2) AXX - The user's box number 
(3) FILENAME - The name of the user's input data deck 
(4) XXXXX - The name of the restart dump file. The user may specify 


any five character name. 
{5) J- HSP sentinel to determine HSP disposition. 
K - Keep HSP files under user's I. D. number 
H - Send HSP files to NIPS printer. 


The transport code will automatically edit to TTY every 1/20 of TMAX (total 
simulation time) unless edit times are specified in the input. The variables 
edited are DELT, TIME, and NTI ‘number of time increments). 


RUNNING CODE FROM ORDER BATCH PROCESSING SYSTEM 


If the user wishes to run the MATTRAN code under the ORDER batch processing 
system, the following control] cards are required: 


(1) ¥*1D ACCOUNT NUMBER NAME BOX AXX 
(2) *ROFILES HL524 VELOUT 

(3) *ADT *AANNN TPPXXXXX TRMXXXXX XXXXX TEPXXXXX TSFXXXXX 

(4) *XEQ MATTRAN 

(5) *DATA 

(6) *AXX FILENAME 

(7) XXXXX J 

(8) special instruction card; blank if not a restart Format (4F 12.4, 15) 


In the above example, tape HL524 contains a file required by the transport 
code; in this case, it is file VELOUT. This file contains the Darcy 
velocities generated by the TAYLOR ground-water code(GWFLC#). Tape AANNN 
line(3) is the user‘s private tape. Files generated by the transport code 
will be written to this tape. The * before the tape name instructs the 
computer operator to hang a blank tape and label it accordingly. Under the 
ORDER system mode, it is imperative that the HSP sentinel on Vine (7) be set 
to H and not to K, Sentinel H will send HSP files to the NIPS printer. 
Option K will keep the files under the I. 0. number and this may result in 
loss of HSP output if the job is being run by computer operations. Special 
instruction card, line (8), contains provisions for changing the values of 
DELT, CHNG, DELMAX, TMAX, and NTI on restarts. They are entered in the above 
order; enter a zero for those values that do not have to be changed. 

The preceding set of control cards are added to the front of the MATTRAN data 
deck and submitted to production. Contact Miriam Lohmann for further details 
on K-Division production job procedure. 
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FILES GENERATED BY TRANSPORT CODE 


The *ADT order control card (line 3) will write the following files to 


magnetic tape. 


f1) 


(2) 


TPPXXXXX 


TRMXXXXX 


XXXXX - 
TF PXXXXX 


TSF XXXXX 


A short description of each file is added. 


This file contains concentration values generated by the 
MATTRAN code. This file is tater postprocessed and 
contour plots are generated. 

This file contains specia! messages from the transport 
code. If running from ORDER, a problem edit 1/20 of TMAX 
will be included unless edit times are specified in the 
input. 

Restart dump fle generated by code. 

Contains variables FRATE(3), FLOW(3), TFLOW(3), FRATE(6), 
FLOW(6), TFLOW(6) and TIME at each iteration. These 
variables can be plotted against TIME by the system plot 
routine PLOTPK upon completion of the run. 

Contains variable arrays generated jin the subroutine 
SURF. If in a preceding run the boundary conditions are 
identical to the problem to be run, this file can be read 
by MATTRAN to eliminate access to subroutine SURF. Ona 
very large problem, as much as 3 minutes may be saved. 


li 


Distance -- m 


EXAMPLE PROBLEM EXECUTED FROM A REMOTE TERMINAL 


To complete this run of the transport code, 6 files are needed: 1) MATTRAN, 
the executable code, (2) UXTV, £3) UXRJ, (4) TRANPP, the gostprocessing code, 
(5) VELOUT, the velocity input file and, (6) the data file, MATTRAN, UXTV,> 
UXRJ? and TRANPP are public files. VELOUT is a velocity file created by 
owFLow’ +6 and functions as an input file to MATTRAN. Because of this, an 
identical mash description must be used as input to both GWFLOW and MATTRAN. 
VELOUT is a binary file that is compatible with MATTRAN. An example problem 
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taken from a material transport study” is presented to demonstrate the 


creation of the data file. 


Figure 1 illustrates the finite element mesh of the example problem, Node 
numbering proceeds along vertical mesh lines from bottom left to top right. 
Element numbering is similar. The problem configuration contained 1562 
elements whose shapes were prescribed by the coordinates of the 1656 nodes 
that formed the corners. 


1 Seepage 


Distance — km 


FIG. 1, Finite element mesh of example problem. 
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The data file of the problem is found in Appendix A. The first line contains 
the problem identification number and title. The second line is the controi 
line. The control line indicates 1656 nodes, 150 elements, 22 different 
materials, a lack of elements with corrected material properties, 10,000 as 
the maximun number of time increments, one Dirichlet boundary condition (the 
source node}, +76 element sides on which flux-type houndary conditions are 
app’ ied !zero flux boundaries’, 10 elements having seepage surfaces {surfaces 
throuqdh which material is released from the system). Also, velncity input 
cantro? equals one ‘which means read the VELOUT file’, storage contro! equals 
zerg ‘this cantrat his hean retained far possible future use and the code has 
heoan aqdif ied to automatically store the output on tisk), KSS equals one ‘to 
ohtain 2 teansient ca*ytian!, and ISSENT is set to 7 ‘in order fnr the SURF 


sukeoutine to caleulate TSFAREXKY. 


Line 3% ig the sinmylatton time cantral line, The control parameter for saving 
roncenteations for the postorocesser “TRANPP’ is included as the first varia- 
hte in Tine 3.) Tt was sot at 9 to produce concentration contour plots at the 
times ‘yeacs) entered in tines 4 and 5. The second variable in line 3, the 
time increment ‘DELT}, 15 9.2 years; the multiplier for increasing DEL’ is 0.1 
years; the maximum yviiue of DELT is 55 years; the maximum simulation time 
‘TMAX' af the prohlem is 1000 vears; and the time integration parameter is 
one, Lines 4 and 5 list the years at which concentration values are saved on 
dist. Line 6 contains output control information, The number 4 ‘s a flag 
commanding the initial corditions t be printed. The number 2? commands fluxes 
ant concentrations to be orinted at the simulation times listed in lines 4 and 
‘ 


y 


. ‘Line 7? indicates one Dirichlet boundary node and line 8 indicates tha: 
this boundary will be chanjed to a new value at 100 years. Lines 9 through 52 
define the properties of the 22 materials, Although .:e capability of 22 


materials exist, these data describe a 5-layer system.*. 


“or ayampte, lines 9 and 10 describe material Tt: the distribution coefficient 
‘quals O, the bulk density is 1.75, the longitudinal dispersivity is 50 m, the 
“ransverse dispersivity is 5m, there is no radioactive decay, porosity equals 
1, the modified coefficient of compressibility is 2ero, moTecular 


diffusivity is lm, and (on line 10) tortuosity is 3.1. Lines 53 through 32% 
describe the noda! position and Dirichlet concentration-type boundary 
conditions. Some of these data are formated and some are comma delimited 
because data generators and text editors were used when convenient. For 
example, line 53 describes node | by: node number equals J and a Tack of 
Dirichlet boundary generation. The second 0 means this is not a Dirichlet 
boundary node (1 means it is a boundary node, e.g,, node 586}, the x 
coordinate position is 0.0 m, the y coordinate position is 1000 m (Fig. 1), 
the Dirichlet boundary value is 0.0 (the concentration at the source node 536 
is 1000.0). 


‘ines 322 through 1883 define the elements. Line 322 describes element one 
as: element number equals 1, the lower left node of the element is 1, the 
lower right node is 24, the upper right node is 25, the upper left node is 2, 
the element properties are those of materia! 1. Lines 1884 and 1885 set the 
initial concentration values at all nodes equal to zero. Lines 1886 through 
2061 define the Neumann flux-type boundaries. The example problem (Fig. 1) 
was designed to release material only at the downgredient end via seepage 
surfaces. Therefore, zero flux-type boundaries surround the remainder of the 
model, These boundaries (element sides) are defined hy two nodal points. One 
such houndary is defined in line 1886 as being between nades | and 2. A zero 
follows (indicating a lack of automatic generation); the final zero in the 
Vine defines the outward flux as zero. The seepage surfaces (element sides) 
are designated in lines 2062 through 2071. In line 2062 the element 1562 has 
a seepaae surface between nodes 1556 and 1633. The zero indicates lack of 
automatic generation. 


After execution, the high speed files (the printed output) contained the input 
data from the data file, VELQUT. This is followed by the material flux 
concentrations and a systems flow table printed for the years in lines 4 and 
5. The graphics output file {in this prob'em called TPPTG231) was used to 
produce the concentration contour plots in Appendix 8. TRANPP was the 
postprocessor used to produce the plots in Appendix B. 
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TRANPP--A POST PROCESSOR FOR THE MATTRAN CODE 


INTRODUCTION 


TRANPP is a LRLTRAN code written to generate graphics from the concentration 
file written by the MATTRAN main code. The TRANPP code uses system libraries 
ORDERLIB, TV8OLIB, and UX80LIB to generate contour plots. A modified version 
of ORDERLIB routine, CONTUR, is used for contouring. The MATTRAN main cade 

weites ta disk at the specified simulation times. The TRANPP postprocessing 


code reads the sets of concentrition values and generates the plots. 
HOW TO RUN THE TRANPP CODE 
The following files are needed to generate contour plots: 
(1) TRANPP - Executable control ‘ee 
(2) TPPXXXXX - A file written by the MATTRAN code containing 


concentration values and other required variables. 


One of the following UXDD80 system files is required and the other two are 


optional: 

(1) UXTV - A routine that is a member of the UXDD80 graphics system; 
jt ts used to view plots on TMDS screens. 

(2) UXRd - Also a member of the UXDD80 graphic system. Tnis routine 
is used to send the UXDD80 plot file to the RJET printer. 

3) UXFR - A UXBD80 system file that can be used to send the file to 


other output devices. 
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I.0. onto a teletype and proceed in the fo‘lowing manner (the following 
parameter values were used to produce Appendix & that contains the graphics of 


the example problem): 


USER: TRANPP / 3 .39 
TTY: TYPE CONCFIL NAME 
USER; TPPTGB31 
TTY: TYPE-NLEV,GRID OPT, XMIN,XMAX, ZMIN, ZMAX, CONC-MIN,CONC-MAX... 
FORMAT(215,6F 10.4) 
USER: 11,0,27490. , 31400. ,0., 1000. ,1.9,1009. 
TTY: DISK FILE NAME. 2 2 wee TPPXXXXX 
# OF CONTOUR LEVELS 2... 2 + ee 11 
GRID OPTION 2. 1 2 1 ee ee ew 0 
AMIN oie 5. care ioe. fon Sse Nae Aone el oe 2.7490E+04 
AMAR. oy a a. Ga ae oe. Se Sees 28 3.1400E+04 
ZEMIN a es oe ee ee ew 0. 
EMAXK cee: Be 5) Bok ee A a Oe 1, 0000E+03 
MIN CONCENTRATION . 1... . . » 1.9000E+00 


MAX CONCENTRATION . . . «~~~ © » 1,0000E+03 


TTY: IF ALL TTY INPUT OK TYPE Y OR HIT RETURN TO REPEAT 
USER: Y 

TTY: TYPE FIRST SET OF LEVEL NUMBERS..10F8.4 
USER: 900. ,500.,100.,50. ,20.,10.,8.,6.,4.,3. 

TTY: TYPE SECOND SET OF LEVEL NUMBERS. .10F8.4 
USER: 2.5 

TTY: THE 20 LEVELS OF CONC VALUES ARE 


900.00 500.00 100.00 50.00 20.00 10.00 8.00 
6.00 4.00 3.00 


2.00 -0. -0. -0. -0. -0. -0. 
-0. -0, -0, 

TTY: ALL GOOD DATA PLOTTED....PROBLEM DONE 

TTY: ALL DONE 
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Having completed execution of the postprocessing controllee, TRANPP, the user 
“ill fave a file or family of files starting with the name UATV80DD80. The 
plot file (UATV80DD80) contains contour plots. At each specified time, the 
code will generate a contour plot with nodes, and a contour plot without 
nodes. A maximum of 40 contour plots are possible. The user may now view the 
plots in file UATV800N86 on TMDS with routine UXTV or send the file to the 
desired output device with routines UXRJ or UXFR. Refer to Blair's report? 


for more information about UXTV, UXJR, and UXFR. 
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CONCFIL ~- 


NLEV ~ 

GRID OPTION - 
0- 
1 - 

AMIN, XMAX - 


ZMIN,ZMAX - 


CONC-MIN - 


CONC-MAX - 


OEF INIT EC OF TERMS 


A file written by the MATTRAN code. It contains 
concentration values plus other variables needed to 
generate contour plots. The file will be of the form 
TPPXXXXX. 

Number of contour levels. This is an integer greater 
than or equal to one, and less than or equal te 20. 
An integer with value of 0 or 1. 

This option wili generate a i page plot. 

This option will generate a plot which is dimensioned 
one-to-one in X-direction vs Z-direction. 

The user's minimum and maximum Cartesian coordinates, 
X-axis. 

The user's minimum and maximum Cartesian coordinates, 
Z-axis. 

The minimum concentration value the user wishes to 
contour. All values less than this number will be set 
tozero internally by the TRANPP code. 

The maximum concentration value the user wishes to 
contour. All concentration values greater than this 
value will be set to zero internally by the code. 
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APPENDIX A: 
INPUT FOR THE EXAMPLE PROBLEM 


This appendix contains data used as input for the example problem described in 
the text of tne attached report. 
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Wr 


OL. INPU PILE FOR EXAMPLE IN TRANSPORT CODT USER“ HAMRAL 
1056. 1562.27-0-10000.1.178.10.1.0.t. 

PaO BID. wid. 510. 

B30.-770, .900. . 100. .0. 


\ 

0. 4295.50.05. 00.0.8. 
gt. 75.50. .5..0...1.0-.8. 
ort. 75.50. .5..0...1.0..1, 
a3. 75.50. .5..0...1.0..0. 


Ot 75.59. 5, .0,..0.0..0. 


0114.75.30. 5.00. ..08.0..1. 
0111.75.30..5..0...05.0..t. 
011.75.30..5.-9...08.0..1. 
OS 1.75.30, 05.0. 05.0.8. 
0151.75.30. -10..0...01.0..1. 
3075.50. .10..0. 010.04. 
0151.75.30. .t0.-0...01.0..1. 
Bt.75. 30.00. 0, Ott. 
D1 1y.75.50..10..0...01.0..1. 
0111.95.30. 5.0. .-05.0..1. 
Bo 11.75.30..5..0...05.0..1. 
0. 69.75.30..5-.0...05.0. 01, 


0..1.75.30,.5.,.0-..05.0..1. 


1 
3..1,.75.50,..5..3. .0.0..1. 
Wd 
0..1.75.50..5..9. .0.0..1. 


1,1 
0,.1.75.50..5.-0-.-1.0..t. 
wt 


0,.1.75.50..5..0...9.0..0, 
ima 

10 Q0,.100D. .0. 

20 00,.101..0. 

23.0 Of..1308..0. 

a4 Oo 

P60 

46 0 

47 9 

48 0 
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69 0 O00O..t188..0. 
79°90 OQF090..820..0. 
71-9 06000. .834..0. 
9@ 9 06000..1128,.0. 
930 OB0td..760..¢. 
oH 9 06000..775..0. 
16 9 QB000..1100..0. 
''6 O O9900..730..0. 
117 9 D9000..746..9. 
'26.9.0.9000..990..0. 
127.9.0.9000. .906..0. 
128.9.0.9099. .972,.9. 
138.9 09000..t082..0. 
13@¢ 9 g10000..700..9. 
wus 9 OlOGNG. .717. .e. 
1619 10000. .1074. .0. 
162 OQ 011000..679. ,0. 
1683 0 011009. .686..9. 
154.0.0.11909..706..0. 
155.0,0,%1009.,7e4, .0, 
165.0.0.t1900. Jer. .0. 
ay 0 911000... °055..0. 
145 0 017900. .640..0. 
1eE 9 012900. .659. .0. 
187.¢.9. 17000. ,.678..0. 
188.9.9.12099. .657..0. 
169.9.0.12000..716..0. 
e079 Gie008. .t058..0. 
POB O 013900. .610. .0. 
P09 9 013009..630. .0. 
239 0D Of3000..1050, .0. 
23! O 014500. .5a0. .0. 
P32 0 014000..601..0. 
253.0 914006..1042..0. 
@54 Oo 019000. .550..n. 
255 @ 916000. .972..0. 
P76 0 015000. .1034..0. 
277 @ 016000, .520,.0. 
278.0 916000. 542. .d. 
299 0 d16000 00+. .0 
300 0 917000 
30' Go a17000 
3ee 2 917000 
3°37 0 O1a0cg. .460. .0. 
Fee O 016000. .4a4. .c. 
345 0 016000. .9aa. .0. 
346 0 Ot9G00..430..0. 
347 9 4819000. .455..0. 
368 0 0190090..9B0..0. 
369 0 Of0000. .400..0. 
370 0 OFOG09, 425. .0- 
380.0.0.20000..666. .0. 
361.0.9.20000..712..0. 
3o°.9,0,°0900. .739..0. 
391 9 Of0000..972..0. 
392? 0 021000. .370..0. 
393 0 O°1000..397..0. 
414 0 OFt000. .S6¥..0. 
415 9 O7000..340..0. 
416 9 Jeenod. . 368. .0. 


“37 0 O22000. .936..0. 
4380 OP3000. .3t0..0. 
433 0 OPxd00. .339..0. 
460 0 023000, .948..0. 
4610 024000, .7B0..0. 
46? 0 pho ee 
4830 Ofxg00, .940 
“8+ 0 029009. .250 
#85 0 Ge5000..281 
S06 0 029000..932 
507 0 Op6a00, .220 
s0a 0 GeB000. .282. 
Se9 0 476000. .9¢4 ‘i 
530 0 Qo. 
Sst o 0. 
552 0 G27000..916..0. 
553 0 0e7800..175..0. 
554 0 027509, .208. .0. 


561.0.0,27500. 929. .0. 
S6?.0,0.27500. .972. ,0. 
563.0.0,07500..505..0. 
564.0,0,27500. $30. .9. 


375 0 De7900..901..0. 

576 0 077750. . 167... 

577 0 027750. .200. .0. 
585.0.0.27750. .464..0. 
1.27750. .497.. 1000, 
0.27750. .53! 
0.27750. . 728. .0. 
0.27790. .751..0. 
Q 
cot 


227750. . 794. 0. 
027750, .893. .f, 


6to.c, '0.26000..523..0. 
6t6.0.0.28000..72t. .0. 
617.0.0.28000. . 75+. .0, 
618.0. +: 28000. .767..0. 


6e1 078000. .686..0. 
672 H Oresog..U52..0. 
673 0 OPasOD..185..0. 
644 0 O7aso0. .B78..0. 
645 0 O2g00D. .tN5..0. 
646 0 OPgN0o. .t7B. .0. 
657 0 079000. .071..0. 
BBB 0 030000..135,.0. 
569 0 030000..169..0, 
590 O O30000..06%..0. 
591 0 O31900..129..0. 
Bye OG 31000..1Ge..0. 
AZO 031000. .655. .0, 
714 9 DBz000..100..0. 
715 0 OSP000..t34..0, 
736 9 O3P7000..a48..0. 
737 Q 133000, .97. .a. 


138 
769 
765 
76! 
7F 
mith 7 
18+ 
eo5 
a06 
B07 
aPb 
ara 
a39 


9 
6 
9 
o 
13) 
(2) 
1) 
ti) 
9 
Q 
9 
0 
9% 


O3IOOU. 17. .0. 
O73000. .845..G. 
DIVONOD. .95. .0. 
O34G00. .129..9. 
D3~O00. .G43..0. 
035000. .95..0. 
935000. .1P9 .D. 
O35000. .843..0. 
035500. 95. 9. 
935500, .129..0. 
B35500. 843. .0. 
035750. .97. .9. 
U3SSTSO, 196. .D. 


838.9.0.35750..398 .9. 
839 .9.0. 35750. .w3e. 0. 
840.9.0. 35750. 466. .0. 
441 .9.0. 359750. 500. .0. 


5) 

Bee 
653 
orm 
875 
e716 
ao7 
e938 
a99 
929 
ge) 
922 
gus 
gus 
545 
966 
967 
368 


gq 
9g 
ss] 
9 
g 
i) 
9 
g 
9 
Q 
9 
9 
Uy) 
) 
9 
Q 
0 
9 


035757. 660. .0. 
936009, .92..0. 
936090. .126..9. 
Q3G099. Gv. .9. 
Q36590..90..0. 
O3GSOD. 104. 0. 
036500. .636..0. 
O37000. .90..0. 
037000, ,104, 9. 
G357000. G38. .9, 
038000. .87..0. 
058080. ,3°3..9. 
030000. .835. .0. 
939000. .85..0. 
D39D00. .119. 0, 
9359000. .833. .. 
941000. .80..9. 
D43900, 54.9. 


979.0.0.41000..'Be..0. 
9721 .9.9.41000. 216. .0. 
972 .0.0.41000. 259. ,0. 
986.0.0.41008.,726..9, 
987.0,0.41000..760..0. 
988.9.0.41000..794. ,0 


989 
ggo 
99! 
Wwe 
10s 
guy 
1935 
G36 
1937 


1946.9.0, 
9ug.0.0. 


a 
9G 
i) 
3 
° 
9 
Qo 
Q 


Q 


D¥1009. G28. .0. 
G43006. .35. 
043000. .#09..D. 
Qu 3900. .623..B. 
gusood..70..0. 

Q45000.. 194. .9. 
045000. .818..0, 
G4u7008. .65..0. 

047000, .99..9. 

#7000. .375..0. 
67900. .405..0. 


1947.9.9.47009. .439..9. 


1958 
1959 
1969 
19Bt 

1908? 
1963 
1198 


epeeescoo 


047909. .813..0. 
949900. .66..5. 
Q49009. .94..0- 
949090. .G08. .0. 
049995. .59..0. 
949995. .93. .9. 
349995. .897..0. 
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1105 0 

1106 6 

1127 6 

1128 0 

1129 6 

1150 0 is 

V151 0 053000. .95. .9. 
(162 D 053000. .89. .0. 
tt73 0 053000. .@03..0, 
(194 O =095000..5%..0, 
(175 Q 095000. .05..0, 
1196 0D 055C000..799. .0. 
(197) OS99TS.. .G. 
1198 0 055975. .84..0 
1219 0 e° «0, 
1220 0 OS56025.. . 
1221 0 058025. .04..0 
12 0 058025. .798..0 
1243.09 057009. .49..0 
tee+ O 057000. .63..0 
1265 9 657000. .797..0. 
(266 0 059000. .45..0. 
1267 O 059000. .79..0. 
1208 0 039000. .793..0. 
1789 0 OG1000. .40..0. 
1290 O 061000. .74..0. 
1311 0 G61000. .768. .0. 
1312 O 063000. .35. .0. 
1313 0 063000. .69..0. 
1321.0.0.6300D. . 341. .0. 
13a2.0.0.63000..375. .0. 


1323.0.0.63000. .409. .0. 
1334 0 063000. .783..0, 


1300 0 057000. .773. 
1381 0 089000. .22.. 
1382 9 069000. -5B. 


wos 0 

1404 O 

1405 0 

1s2e6 0 

we? O 

12a 0 

1449 0 . . 
1450 G 075900. .19..0. 
1451 © O75909. .44. .0. 
wre Q O75000..758..0 
1473 @ 077900, .9..0. 
1474 @ O77G00. .42..0. 
1495 0 O77000..756..0. 
1496 0 O77500..5..9. 
W497 0 OT7S00. .39..0. 


1505.0.9.77500..311..0. 
1596.0.9.77506. . 345. .0. 
1507.0.0.77500. .379..0. 
vie 0 ITTSOC. .753..0. 
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1514 0 O7B000..4. 
15PG 0 978000. .38..0. 
641 0 076000. .7527. .0. 


‘42 0 078500..2.. 


1543 G O78500. .36 


4559.0.9. 78500. 689. . 
1569.0.9, 78500. .6%4.. 
1561.0.0.78500. .646. . 


1564 
1965 
1566 
1997 
1966 
1989 
1610 
1611 

1612 
1634 
16 
1635 
1656 


eervFescoocee 


O78500, . 7%! 
O79000..1. 


979000. .35. .0 
Q79000. .749..0. 
o79500..0..9. 
079500. : 
079500. 
Q7978D. .0..0. 
079750. . 
079757. . 748. .0 
geoeo). .2..0. 
OB0900..34..0. 
080000. . 748. .0. 
hs e4 a5 e 
2 PS 26 3 
3 26 eT + 
u 2? 2 5 
5 78 2 6 
6 2; 30 7 
7 30 ci! a 
a cil cd 3 
a Be 33 ig 
10 FF Ww 
vw i bd Ve 
we 3 3 13 
{3 36 0 4 
4% 37 88 15 
19 7a 39 16 
16 39 40 IT 
17 4 1 1B 
1B 41 4 19 
1% we 4 eo 
20 43 44 et 
Cal 44 45 Be 
ee 45 46 23 
e4 47 «48 (PS 
r-] 48 4Q ?6 
e6 49 «6450 ay 
27? 650)06=«(45l 28 
238 BI 5e 79 
29 Se 4463 30 
30 53 5 EM 
Fi) 54 ii] xe 
ye 86665 5S 33 
B 56 57 Fi] 
34 57 58 35 
ci] 5a 59 36 
36 59 «(#6 ” 
Ey 69 61 38 
3 #61 62 39 
33 Ge 6348 
4“) 63 64 “1 


ee ee ee ee ee ee ee 


eonpVvccCoCo CCC CVnoOFDOUOeOKSaD 


o 


0. 
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1907 
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1410 
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116 
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Wg 
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lee 
124 
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139 
wy 
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a7 
149 
14D 
re) 

ee 
143 
144 
145 
146 
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APPENDIX B: 
GRAPHICS FROM THE EXAMPLE PROBLEM 


This appendix contains the graphics generated by the example problem described 
in the text of the attached report. 
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